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Solution of the Bohr hamiltonian for soft triaxial nuclei 



L. Fortunate*, S. Dc Baerdemacker, K. Heyde 
Vakgroep subatomaire en stralingsfysica, Proeftuinstraat 86, B-9000, Ghent (Belgium) 

The Bohr-Mottelson model is solved for a generic soft triaxial nucleus, separating the Bohr 
hamiltonian exactly and using a number of different model-potentials: a displaced harmonic os- 
cillator in 7, which is solved with an approximated algebraic technique, and Coulomb/Kratzer, 
harmonic/Davidson and infinite square well potentials in f3, which are solved exactly. In each case 
we derive analytic expressions for the eigenenergies which are then used to calculate energy spectra. 

Here we study the chain of osmium isotopes and we compare our results with experimental 
information and previous calculations. 
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PACS numbers: 21.10.Rc,21.60.Ev 



I. INTRODUCTION 



The Bohr-Mottelson collective model of the nucleus 
has recently attracted a significant interest because of 
the possibility to derive many "new" solvable cases [1- 
18]. Intense efforts did arise because of the availability 
of models, based on the square well potential, that are 
related with the issue of critical point symmetries at the 
shape phase transition (E(5), X(5) and Y(5)) [1-3]. This 
has given rise from one side to many applications aimed 
at the survey of existing experimental spectroscopic data 
and at the identification of signatures for the new mod- 
els [19-25] in various mass regions, especially in connec- 
tion with the effort to build new descriptions for tran- 
sitional nuclei. On the other side a number of mathe- 
matical solutions of the Schrodinger equation associated 
with the Bohr hamiltonian with various model-potentials 
have been proposed [5-11, 14-18]. For some of those po- 
tentials, that in general are functions of the quadrupolc 
deformation variables (3 and 7, an exact separation of 
variables is possible, while in other cases an approximate 
separation holds. Usually the spectrum (and transition 
rates) is derived analytically and compared with available 
experimental data. We like to mention that for a number 
of these new solvable cases extensive comparisons with 
experimental data are still not available. 

Recently a solvable model was proposed for the soft 
triaxial rotor with a minimum in the potential along the 
7 direction located at 7r/6 [14]. The aim of the present 
paper is to extend and complete that solution, propos- 
ing a solvable model for which the Schrodinger equation 
is separable. The 7— angular part is then solved for a 
harmonic potential centered around any 70 in the inter- 
val ]0,7r/3[ , while the (3 part is solved exactly for the 
Kratzer-likc potential, for the Davidson potential and for 
the infinite square well potential. 
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FIG. 1: Contour plot of the potential V(/3, 7) with minimum 
in 70 = 10° and flo = 0.2. The situation depicted here cor- 
responds to a very soft deformed triaxial rotor, where the 
Vi(f3) potential (see Eq. (3)) has a Kratzer-like form, while 
the Vi (7) is a displaced harmonic oscillator. The coordinates 
are X = /3 cos 7 and Y — /3shi7. In the point /3 = the 
potential goes to +00 (black area). 



Another reason for the search of analytic solutions for 
the more general case with non-irrotational moments of 
inertia comes also from a recent study [26] in which the 
rigid triaxial rotor model is improved by relaxing the as- 
sumption of irrotational flow moments of inertia. The 
fit of the three components of the moment of inertia to 
experimental data improves the description of a number 
of nuclear properties and suggests that the irrotational 
assumption is not correct. Our model incorporates this 
idea and extends the rigid model to 7-soft models. 

To give an idea of the typical problem that we are con- 
sidering we show as an illustration in Fig. (1) a contour 
plot of the potential surface with a Kratzer-likc potcn- 
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tial in (3 (with the minimum at (3q = 0.2) and a dis- 
placed harmonic oscillator potential in 7 (with the mini- 
mum at 70 = 10°). The solution of the two-dimensional 
Schrodinger equation for the Bohr hamiltonian contain- 
ing this potential surface is given by a set of eigenval- 
ues and their corresponding cigenstates. To each eigen- 
state (and in particular to the ground state) an average 
quadrupole shape is associated. In the present case (Fig. 
1) the wave functions are concentrated around the re- 
gion of the minimum and thus we are dealing with a soft 
triaxial rotor. 

The present study may also be considered as an at- 
tempt to solve a simple model that may be compared 
with microscopic theories (like HF, HF+BCS or HFB) 
[27-29]. The potential surfaces used throughout the pa- 
per, rather than calculated with a variational procedure 
in a microscopic framework, are instead a priori given 
in analytical form. A number of early microscopic stud- 
ies, like for example Ref. [27], concluded that 7— softness 
was not only unavoidable in triaxial nuclei, but that the 
potential barrier in the 7 direction was so small to cast 
doubts on the applicability of the rigid rotor model. We 
therefore reanalyze a number of osmium isotopes, which 
were treated within a HFB formalism, and compare our 
results with the results of Ref. [29] . 

In the present paper, we use a potential of the form 
Vi(j3) + V2(7)//3 2 , which allows for a separation of the 
Bohr hamiltonian in (3 and (7,6^) part (Section II). The 
solution of the (7,^) part is amply presented in Section 
III. In Section IV we discuss the quantum numbers asso- 
ciated with the present model and presents a few com- 
ments on quasidynamical symmetry. In Section V we dis- 
cuss in detail some examples that illustrate the method 



of solving Eq. (10). In Section VI, we solve analytically 
the (3 part of the problem in the Coulomb/Kratzer, har- 
monic/Davidson and infinite square well cases. In Sec- 
tions VII and VIII, we discuss the fitting procedure and 
its accuracy and we apply the method to the chain of 
osmium isotopes and make a comparison with other cal- 
culations. Here, we also stress the fact that, in doing 
so, in general we go beyond an irrotational approach and 
make use of three different moments of inertia used as 
parameters. The need to take into account this possi- 
bility has been emphasized before by Wood et al. [26]. 
We also compare with the situation of irrotational flow 
in which a value 70 can be deduced and implies that 
fluctuations of the moment of inertia in the 7 direction 
are neglected, but the softness is included. It is interest- 
ing to notice that the values of 70 so obtained correlate 
well with results extracted from different methods used 
to analyze data in the Os nuclei (Section VIII). We ana- 
lyze these isotopes, particularly because they are thought 
to be situated in the transitional region between 7— rigid 
and 7— soft shapes. In Section IX, we formulate the main 
conclusions of the present study, while in appendix A, we 
explicitly present the matrix elements used in the calcu- 
lations. 

II. FORMULATION OF THE MODEL 

The Schrodinger equation for the Bohr hamiltonian 
reads 

ff fl tf(/3, 7 A) = £*(/3,7A), (1) 
where the hamiltonian is given by 
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Here Li are the projections of the angular momentum (L) 
on the body-fixed axes and B m is the mass parameter. 
An exact separation of the variables (3 and 7 may be 
achieved when the potential is chosen as 

VG9.7) = + (3) 

The resulting set of differential equations (one contain- 
ing only the (3 variable and the other containing the 7 
variable and the Euler angles, L with l = 1,2,3), after 
multiplication by 2B m /h 2 , reads 



I sin (3 7 )^ Sm( ^7 + 2^1 3 4 sin2 (7 ~ 27 ™/3) 

+ua(7)-w}#y,0 t ) = 0, (5) 

where co is the separation constant, e = (2B m /h 2 )E and 
u t = (2B m /H 2 )Vt with t = 1,2. The wavefunctions sat- 
isfy #(/?, 7, 6 L ) = £(/?)V>(7A)- In Ref. [14] the same 
derivation has been proposed, but the problem was re- 
stricted to deriving the solution of the soft triaxial rotor 
around 7 ~ 7r/6. In that case the rotational kinetic term 
may be written in a very simple way and the solution of 
7— angular part may be given in a straightforward way. 
The main result of Ref. [14] was, (i) to extend the Meyer- 
ter-Vehn formula [30] (strictly valid for a rigid rotor at 
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7 = 30°) proving that the corresponding formula for the 
7— soft rotor requires the addition of a (trivial) harmonic 
term in the 7 degree of freedom and (ii) to show that the 
solution of the full problem does not necessarily imply 
a trivial extension with some additive term to include 
the (3— vibrations, but rather yields an expression for the 
energy levels in which the quantum numbers are inter- 
twined in a more complicated way. 

Here we extend the results of Rcf. [14], addressing the 
more general problem of a soft triaxial rotor (not con- 
fined to 7 ~ 30°), solving approximatively the 7— angular 
equation with a potential of the form 

w 2 (7)-C'(7-7o) 2 = Cx 2 , (6) 

that has a minimum in the interval 0° < 70 < 60°. For 
symmetry reasons we can restrict ourselves to the sector 
0° < 70 < 30°. The obstacle to the further separation 
of variables in Eq. (5) is represented by the rotational 
term (second term), that mixes the variable 7 with the 
projections of the angular momentum. This term can be 
rewritten as J^iAiLj. We will follow two strategies to 
deal with the coefficients Ai. 

• The first strategy is to approximate these coeffi- 
cient as follows 

At = 5 , (7) 

4sin 2 (7o - 2iri/3) 

and to replace the corresponding terms in Eq. (5). 
The physical meaning of this approximation is that 
the fluctuations of the moments of inertia are com- 
pletely neglected, while the softness is taken into 
account. Making use of the form (7), we can 
parametrise the components of the moment of in- 
ertia by one single parameter. 

• It has been shown in an extension of the rigid rotor 
model [26] that the components of the moment of 



inertia fitted to experimental data do not necessar- 
ily agree with the irrotational approximation. We 
will therefore include this observation in our model, 
as a second strategy, taking the components of the 
moment of inertia as parameters. 

In the spirit of finding a simple solution, we introduce 
in Eq. (5) the further simplification sin 37 ~ sin 370 , 
obtaining 

{--^+Cx 2 -u + AiL*}i>(x,9)=0, (8) 

i=l,2,3 

where the variable x, introduced in (6), has been used 
and where the second order differential operator D 2 has 
been defined. The wavefunction may be written in the 
unsymmetrized form as an expansion in terms of rota- 
tional wavefunctions, namely 

K' 

where K' is introduced to distinguish the multiple occur- 
rence of states with the same L. One may notice that the 
rotational basis has an SO(2) symmetry reflected in the 
good quantum number K' . We now multiply to the left 
by T>\ L k an d integrate over the Euler angles, exploit- 
ing the orthonormality property of the Wigner functions 
where possible. The result is the set of equations (one 
for each allowed value of K) 

(D 2 -cj)a 1 j' { {x)+J2Aj2( LMK I £ i I LMK')a L K ,{x) = 0, 

i K' 

(10) 

each of which contains matrix elements of the squared 
components of the angular momentum (see appendix). 
We may rewrite it also as 



Do - co 



Y^M LMK I 4 ? I LMK)\a^(x)- 



^ — ^-(LMK I L\ I LMK - 2)af { _ 2 (x) + ^ — —^(LMK | L 2 _ | LMK + 2)a L K+2 {x) = (11) 



to highlight that only three of the a coefficients appear III. ALGEBRAIC APPROACH 

in every equation. 

A convenient way to treat the 7— angular part is to 
consider the dynamic symmetry associated with the D 2 
operator, that is essentially a harmonic oscillator hamil- 
tonian. The advantage here lies in the fact that the 
spectrum generating algebra allows us to write explicit 
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expressions for the 7— vibrational spectrum. 

The infinitesimal generators of the sp(2,M) Lie algebra 
may be written as 

(12) 

where A: is a constant. Similar approaches have been 
used in Ref. [4, 6]. The operators above close under 
commutation, i.e.: 

[Zi, Z 2 ] = -MZ 3 [Z 3 , Z 2 ) = -2iZ 2 [Za, Z x \ = 2iZ x . 

(13) 

With the linear transformation 

X\ = —(Z1—Z2) X 2 = -Z3 X 3 = -(Z1+Z2) , 

(14) 

one may recognize the standard commutation 
relations of the four isomorphic Lie algebrae 
su(l,l)~so(2,l)~sp(2,l:)~ sl(2,R) 

[Xi,X 2 ] = -iX 3 [X 2 ,X 3 ] = iXi [X 3 ,Xi] = iX 2 . 

(15) 

It is also very useful to define raising, lowering and weight 
operators for this algebra, 

X±=Xi± iX 2 Xo = X 3 , (16) 

that obey the following commutation relations 

[X+,X-] = -2^o [X ,X ± ] = ±X ± . (17) 

The action of the above operators on orthonormal basis 
states for the irreps of su(l,l) (| n, A) with n = 0, 1, 2, ..) 
is given by the equations: 

X+ I n\) = y/(X + n)(n+ 1) | n+l,X) (18) 
X_ I n + 1, A) = V(A + n)(n + 1) | nA) (19) 

X I nA) = -(A + 2n) I nX) . (20) 

The values of A are found by comparing the standard 
eigenvalue equation for the Casimir operator of su(l,l) 
with the eigenvalue equation for the same operator, but 
explicitly realized in the terms of the operators defined 
in (12): 

C 2 I nX) = (if - X\ - If) I nA) = ^A(A - 2) | nA) . 

(21) 

The action of the Casimir operator on a given function 
(p(x) gives 

C 2 \<t>) = U(^a + Z 2 Z X ) - ^1 I 4,) = ~ I 4>) . 

(22) 

Therefore we obtain A = 1/2,3/2. As a consequence of 
the definitions given in this and in the previous section, 
we notice that 

D 2 = y/CiZi + Zi) = VC4X , (23) 



where the constant k is defined as k = \fC. The action 
of the operator 4X is given, for each allowed value of A, 
by (l + 4n) and (3 + 4n) respectively. Combining the two 
results we may write that the eigenvalue equation for the 
D 2 operator is given by 

D 2 I n 7 ) = VC(1 + 2n 7 ) | n 7 ) , (24) 

where we dropped the index A to cover the whole spec- 
trum, and as D 2 is the harmonic oscillator operator in 
the 7— degree of freedom, we can associate n with the 
number of phonons in the 7— variable, n = n 7 . 



IV. SOLUTION OF COUPLED EQUATIONS, 
QUASIDYNAMIC SYMMETRY, QUANTUM 
NUMBERS AND BAND STRUCTURE 

We discuss here the quantum numbers that appear 
when solving the problem presented in the previous sec- 
tion. The coupled set of equations (10) contains the rota- 
tional and 7— vibrational structure through the functions 
a^-, and the operator D 2 respectively. In the absence of 
7— vibrations (n 7 ~ 0), we are left with the rotational 
structure (see Fig. 2), which can be analyzed as follow- 
ing. For each L there is a set of allowed if's: {-fQ, Kf} 
where Ki = and Kf — L for L even (thus we have 
L/2 + 1 coupled equations) and Ki = 2 and Kf = L — 1 
for L odd (giving (L — l)/2 equations). 

For 7 = and 7 — 7r/3 the projection of the third 
component of the angular momentum on the intrinsic 
axis 3 gives a good quantum number (K), while for 7 = 
7r/6 the eigenvalue of the projection on the intrinsic axis 
1 is a good quantum number (R). In the intermediate 
regions, none of them may be taken as a good quantum 
number. In the expansion (9) the introduction of K' as 
a label to distinguish between the multiple occurrence of 
states with the same value of L is justified because it is 
referring to the rotational states. 

It has been observed that in the triaxial region, mov- 
ing from 7 = 30° toward 7 = 0°, different groups of 
states may well be classified into bands: a first band 
(0 + ,2 + ,4 + , ...) tends to the finite axial rotor values and 
corresponds to K = 0; a second band (2 + ,3 + ,4 + , ...) is 
identified by its behaviour when 7 — > (in Fig. 2 this 
group of states somewhat cluster around 7 ~ 10° — 12°); 
the beginning of a third band (4 + ,5 + , ....) is seen to es- 
cape to infinity at a quicker pace (leaving Fig. 2 at 
around 7 ~ 20°). 

The experimental observation that a classification in (3 
and 7 bands seems an almost universal feature of nuclear 
spectra reinforces this choice. The labeling with the K 
quantum number is often encountered in the literature, 
although for what we have said here it may not be consid- 
ered adequate. We will in the following use the notation 
K* to regroup the various eigenstates into bands, which 
are the counterpart of the bands with good K, that are 
present in the axial cases. 
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FIG. 2: Eigenvalues of Eq. (10) as a function of the asym- 
metry for n 7 = (relative to the lowest eigenvalue). In this 
case the solution of the rotational part for the 7— soft triaxial 
rotor corresponds to the outcome of the Davydov model [31]. 
In addition, the soft model has infinitely many other copies 
(one for each possible value of n 7 ) of these sets of bands at 
higher energies (depending, of course, on the magnitude of 
C). 



structure of the 'maximally' triaxial rotor at 7 = 30° 
persists also in the region around 7 ~ 20° — 30°. In the 
intermediate region these groups of states escape to in- 
finity. It must be further noticed that the region where 
the various states, that originate from the axial rotor 
side, are most affected is exactly the region where the 
states coming from the 7 = 30° triaxial rotor diverge. 
The strange character of this quasidynamical symmetry 
is that (at variance with the cases discussed by Rowc 
and collaborators [17, 32, 33], where a true phase tran- 
sition was present between two exactly solvable limits 
associated with different symmetries and different group 
structures) here we are dealing with a smooth transition 
between two limits which formally have the same un- 
derlying group structure, SO (2), and there is no critical 
point in between. Therefore we conclude that the pro- 
posed labeling, that retains the formal division in /3,7 
and K bands typical of an axial rotor, is not only jus- 
tified by the empirical observation that non-axial nuclei 
display the same classification in bands, but it is also jus- 
tified in view of arguments based on a group theoretical 
approach. It is not clear, at present, if a quantization 
procedure around a tilted axis may help to shed light on 
this aspect. 



V. EXAMPLES 



For phonons with total angular momentum 2, one usu- 
ally takes into consideration two projections | K | = 
0, 2. These projections correspond to the so-called 
/3— and 7— phonons, that generate axial and non-axial 
quadrupole oscillations, respectively. Because of the non- 
axial character of the 7~variable, the 7— vibrations are 
naturally incorporated into Eq. (10). This is obtained 
through the operator D 2 , which adds the 7— vibrations 
to the basic structure of rotations of a rigid triaxial body. 
Consequently, one can construct 7— bands that consist of 
copies of the structure of the rigid triaxial rotor (Fig. 
2), where every copy possesses a different number of 
7— vibrations, characterized by n 7 . Once the rotational 
and vibrational structure is determined, the eigenvalue lu 
can be calculated (in the standard way) and plugged into 
Eq. 4 to solve for the (3— part of the problem. An elabo- 
rate discussion on the (3— part and the chosen potentials 
will be addressed in section VI. 

One may describe the present situation in terms of a 
quasidynamical symmetry [32, 33] of a peculiar charac- 
ter: at 7 = 0° the group SO(2) is a symmetry of the sys- 
tem, associated with K, while at 7 = 30° another SO(2) 
group is a symmetry of the system, associated with i?, 
being the chain U(5)dSO(3) common to the whole sector 
< 7 < 7r/6. In the intermediate region, < 7 < 7r/6, 
the SO (2) symmetry is broken, but it must be noticed 
(see Fig. 2) that the structure of the rotational spec- 
trum, present at 7 = 0°, persists in the whole sector 
without being altered in a dramatic way. Only a smooth 
and slight change may be seen. On the other side the 



Equation (10) is solved here in a few cases. Applying 
the recipe discussed in Sec. IV, the set of differential 
equations is turned in a single algebraic equation in lu. 
When L — 0, only the value K' = is present and there- 
fore Eq. (10) reduces to 

(t> 2 -Lo + ^Ai{QQQ I Lj I 000))ag(a;) = 0, (25) 

i 

where all the matrix elements are evaluated to be zero 
and the first solution is thus WL=o.if*=o,n 7 =o = VC and 
we may write, in general, lol=o,k *=o,n T = \ C{1 + 2n 7 ). 

States with L = 1 are not present in this model. 

For L = 2, the two possibilities are K' = 0, 2, corre- 
sponding to 2 coupled equations, 

(D 2 - w + ZA X + 3A 2 )a% + y/3(Ai - A 2 )a 2 2 = 

(26) 

(D 2 - oj + Ai + A 2 + 4A 3 )al + V3(Ai - A 2 )a 2 = , 

(27) 

whose solutions are, lul=2,k* ,n 7 = 



v / C(l + 2n 7 ) + 2^A 4 ±2 J^A% Y^-\.X- (28) 

i y i i<3 

The rotational parts of these expressions reduce to the 
correct values, 3 and 6 respectively, when 7 = 30°. 
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For L = 3 the only possibility is K' 
tion becomes 



2 and the solu- 



Wi=3,jf *=2,n 7 = + 2n 7 ) + 4A X + 4A 2 + 4A 3 . (29) 

Notice that when 7 = 30° two of the three components 
of the moment of inertia (A2 and A3) are equal to 1 and 
the remaining (A\) is 1/4 so the rotational part of the 
energy becomes 9 (as obtained in refs. [14] and [30]). In 
general we have uJL=3,K*=2,n^ = VC(l + 2n 7 ) + 4 VM !; . 

For L = 4 we can write the determinant of the matrix 
mentioned in Sec. IV as a third degree equation in u>. Its 
three solutions, corresponding to the cases K* = 0, 2, 4, 
may be found analytically, although their expressions are 
rather lengthy. 

It is possible to write an algebraic equation in w for 
every value of L, but this equation may be solved ana- 
lytically only for the lowest values. We must therefore 
resort to numerical computation for high values of the 
angular momentum. 

In Fig. 2 we plot the value of u) for various states for 
C = with irrotational moment of inertia. This case 
correspond to the well-known rigid rotor solution, that is 
a particular case of our model (sec Ch. 9 of Rcf. [35], for 
example) . 

Finally we notice that, as a consequence of the defini- 
tions given above, the following relation holds: 



W200 + ^220 — ^000 + 



(30) 



This expression is a generalization of the well-known re- 
lation E(2+) + E{2' + ) = E(3+) [35]. 



VI. SOLUTION OF THE /3-PART 

Once the 7— part is solved with a particular choice for 
the moments of inertia, we can determine the solution of 
the (3— part inserting the appropriate values for u> in Eq. 
(4). With the substitution £(/3) = x(P)0~ 2 , Eq. (4) is 
simplified to its standard form: 



d 2 xW 
dp 2 



(e-ui(/9) 



(2 + hp 
P 2 



}x(/3) 



= 0. 



(31) 



As shown in [6, 7, 14] an interesting case is represented by 
the Kratzer-like potential (which reduces to a Coulomb- 
type potential when B = 0) 



ui(/3) = -A/(3 + B/j3 2 



(32) 



We depict in Fig. 3, as an illustration, the surface cor- 
responding to the reduced potential u(J3, 7) = u\{f3) + 
U2{j)/P 2 where the potential u\{0) has a Kratzer-like 
form, while the potential U2 (7) has a displaced harmonic 
dependence centered (as an example) around tt/6. The 
depth of the pocket centered around the minimum in 
(/3o, 70) may be adjusted with the parameters of the two 
potentials, A, B and C. 




FIG. 3: (Color online). Plot of tt(/3, 7) as in formula (3) for 
< 7 < 7r/3. The potential m(/3) has a Kratzer-like form, 
while the potential 112(7) has a displaced harmonic depen- 
dence centered around 7r/6 (for purpose of illustration). 



Following [14] we can write the solution to equation 
(4) directly as: 



e(n 7 ,np,L,K*) 



A 2 /A 



^9/4 + B + uj LJ<%rH + 1/2 + np} 

(33) 

where the values of Uh,K*,n have been found analyti- 
cally or numerically as discussed in the previous section. 
Energies are usually redefined fixing the ground state to 
zero and using the energy of the first 2 + state as unit, 
namely 

e(n„n p ,L,K ) = ^ Q> % Q) _ ^ ^ ^ Q) . (34) 

Another interesting case that also leads to a solvable dif- 
ferential equation is the Davidson potential [4, 36, 37]: 



B 



D 



(35) 



This is an extension of the harmonic potential (that can 
be obtained when Bjj = 0). Inserting this potential in 
Eq. (4) gives formally the Laguerre differential equation 
and the spectrum may be written as 



5/2 , (36) 



e D (n^,np,L, K*) = \JAu[2n (} + T L . K , 



where, dropping the indices for simplicity, r is the solu- 
tion of r(r + 3) = Bd + to. The redefinition in reduced 
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energy units takes in this case a very compact form: 
eD{n.y,n/3,L,K ) = . (37) 

72,0,0 — 7"0, 0,0 

The combination of the results for the 7— part in Sec. 
IV-V and the results for the (3— part yields spectra with 
an ample range of different possible behaviours. 

Another interesting case is inspired by the E(5) and 
X(5) solutions [1, 2]. Starting again from Eq. (4) one can 
adopt a different substitution, namely 0(/3) = /3 3 / 2 £(/3), 
and a change of variables, z = These transfor- 

mations, together with the choice of the potential as an 
infinite square well, yield the Bessel differential equation 



0. 



The solution of the equation above is given in terms of 
Bessel J functions of irrational order 

where c s<ul are normalization constants (given analyti- 
cally in [38]). The boundary condition at the wall of the 
potential well, 4>{(3 W ) = 0, implies that the spectrum is 
given by 



esQinp, n~[,L, K*) 



■ 



(40) 



Here and in Eq. (39), x s ,ui is the s-th zero of the Bessel 
function with index that depends on u)L,K*,n^- Notice 
that np = s. The reduced spectrum takes the following 
form 

2 2 

/ r r^*\ X »>UL,K*. ~ ^0,01000 , . , . 

£SQ(np,r>~„L,K ) = - 2 3l_ . (41) 



This solution is somewhat similar to the so-called Z(5) 
solution [11] where an infinite square well in j3 was 
combined with a harmonic oscillator centered around 
7 = 30°. Major differences are the choice of exact separa- 
tion of variables made here, the extension of the solution 
to the whole sector 0° < 7 < 60° and the possibility to 
relax the hypothesis of irrotational motion in deriving for 
the components of the moments of inertia. 



VII. RANDOM WALK FITTING PROCEDURE 

We are now equipped with a procedure to calculate the 
eigenvalues of the 7— part of the problem and we have 
given a few analytic solutions (Coulomb and Kratzer, 
harmonic and Davidson, and infinite square well) for the 
(3— part. The energy spectrum (33) depends formally on 
six parameters, A, B, C, A\ , A2, A3, three of which come 
from the potentials (32) and (6), and the other three are 
the components of the moment of inertia. A closer look 



to Eq. (34) reveals that, once the eigenspectrum is scaled 
in the standard way, it does not really depend on A, while 
(33) does. We can now follow the strategy to keep the 
irrotational hypothesis and, since the various components 
of the moments of inertia are connected to each other by 
means of the relations in (7), one has to deal with just 
one component (or alternatively with 70). The other two 
components may be thus determined inverting Eq. (7): 



7o 



arcsini 

V2VZ3 



(42) 



and, substituting the value of 70 in the definitions of A\ 
and A2, we obtain the results 



(38) Ai = 



4A, 



(y^As 



A, = 



4Aa 



-3 + 1 



V12A 3 - 3 - 1 
(43) 

Therefore our model, using the Kratzer potential, has 
only three parameters, B, C and A3, that may be used 
to fit experimental energy spectra. 

This reduction of parameters, though convenient, may 
be very restrictive (sec [26] ) and we prefer to implement 
also a procedure to keep all the components as inde- 
pendent parameters, at the price of complicating the fit 
to experimental data. We distinguish the fitting pro- 
cedures with the words irrotational (3 parameters) and 
non-irrotational fit (5 parameters). 

Similar considerations apply to the energy spectrum 
of the harmonic/Davidson potential. In Eq. (36) six 
parameters are formally present, but using the reduced 
spectrum and the relations among the components of the 
moment of inertia, one is lead to a dependence on three 
of them only (Br>, C and A3). Alternatively the non- 
irrotational fit includes 5 parameters (Bp, C , A\, Ai and 
A3). Likewise, for the square well case, (4) 2 parameters 
are present in the (non-)irrotational fit. 

Due to the number of parameters we have preferred to 
use a numerical method based on a random walk proce- 
dure in order to determine the energy spectrum. 

As a first step we consider an isotope and a subset 
of its experimental levels (usually quite small, 4 levels 
typically). Starting from some initial set of parameters 
(typically B = C = 20 and A3 = 1), we minimize the 
deviation of the calculated and experimental energy val- 
ues by walking randomly in a suitable part of the many- 
dimensional parameter space. The random walk is ini- 
tially constrained to take only irrotational moment of 
inertia (i.e. A\ and A2 are calculated directly from A3) 
and consists of a coarse-grain and a fine-grain stage, each 
of which takes a certain number of steps. The coarse and 
the fine phases of the procedure differ in the percentual 
amount of the initial parameter values which is allowed to 
change randomly (50% and 10% typically). We usually 
end up with a set of parameters that correspond to some 
local minimum of our hypcrsurface (although it must be 
noticed that in this way we are not sure to catch the abso- 
lute global minimum). From this set one may calculate a 
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FIG. 4: Experimental spectrum of 188 Os (a) and irrotational 
fit (b) obtained with B D = 62.301, C = 82.66, A 3 = 2.3917. 
Experimental data are taken from [39]. 



spectrum, which can be compared with the experimental 



one and a total error, J ^(1 — E\ h / E^ xp ) 2 , that gives 
an idea of the accuracy of the description that one may 
get. Usually, using as input values energies of the ground 
and 7— bands, the typical accuracy may range from 1% 
to 0.1%. At this point one is still allowed to identify a 
value of 70, calculated from the irrotational moments of 
inertia (42), with the angular position of the minimum. 
As a final step we start from the set of parameters re- 
sulting from the irrotational fit, but now we relax the 
irrotational hypothesis, treating all the components of 
the moment of inertia as individual parameters. Often 
the higher freedom of this fit reduces the error by a fac- 
tor 2-10, but sometimes this procedure does not improve 
substantially the quality of the fits. It has to be stressed 
that at this point the identification with a single value 
for 70 does not make sense anymore. 

Our goal is, from one side, to give a simple model that 
goes beyond the estimate of ground state energies and 
furnishes a decent description of ground and 7— bands, 
and, from the other side, allows one to extract a possible 
behaviour of the potential energy surface associated with 
a given isotope. 

One must keep in mind that this model will work only 
when the potential energy surface has just one minimum. 
The fact that, after fitting, some important discrepancies 
result between calculated and experimental values may 
point out that a geometrical model may not be well ap- 
plicable after all. 
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FIG. 5: Experimental spectrum of 192 Os (a) and irrotational 
fit (b) obtained with B D = 132.182, C = 44.6504, A3 = 
1.39954. Experimental data are taken from [39]. 



VIII. THE OSMIUM ISOTOPES 

We applied the procedure explained in the previous 
section to the case of the osmium isotopes with A = 
172 - 180, 184 - 192. The irrotational fit with the David- 
son potential, starting from the 4 + and 6 + levels of the 
ground state band and the 2 + level of the K* = 2, n 7 = 
band and the + level of the K* = 0, n 7 = 1 band, 
leads very quickly to a set of parameters Bd,C,As that 
gives always a rather good description of the low-lying 
positive parity energy spectrum, with an estimate on 
the total error, as defined above, which is in the range 
0.01 — 0.09. In figs. 4,5 we compare the experimental 
energy levels of 188 > 192 0s (a) with the irrotational fit 
(b), obtained from 1600 runs of coarse grained random 
walk plus 1600 runs of fine grained random walk. In this 
case the non-irrotational fit does not improve in a siz- 
able way the quality of the fit. The absolute error in 
the case of the 10 + state is ~ 100 kcV, or less. In the 
best cases the accuracy for the 8 + and 10 + members of 
the ground state band may reach the value 1/1000. It 
is observed that the ground-state and excited bands are 
fairly well reproduced (which is not surprising since we 
used the positions of two levels of each band as an input 
to our calculations). A bit more surprising is the quali- 
tative behaviour of other bands: usually we find that the 
(A'* = 0,n 7 = 1,710 = 0) and (K* = 4,n 7 = 0,71/3 = 0) 
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A Os 

FIG. 6: Reduced energies of the bandheads of the 
(K\ ni ,np) = (2,0,0), (0,1,0), (4,0,0) and (0,0,1) bands 
(indicated with circles, squares, diamonds and triangles re- 
spectively) as obtained from the irrotational fit procedure for 
the whole chain of osmium isotopes. 

bands lie lower (in a few cases dramatically lower) than 
the first /3-band (K* = 0, n 7 = 0, = 1). In Fig. (6) 
we report the reduced energy of the bandheads of the 
lowest bands for the whole chain. In the osmium nuclei 
with lower mass number, we find comparable energies for 
the (3— band and for the K* = band, while the K* = 4 
lies at considerably higher energies. This is reflected in 
the magnitudes of the Bu and C coefficients: for 192 0s 
they are respectively, ~ 130 and ~ 44, but as soon as we 
move to lower mass numbers the value of Bp drops to 
less than 1 for the lightest ones. Thus, the (3— vibrations 
are found at much lower energies. 

The example displayed in Fig. 4 is rather rigid and the 
(dimcnsionless) values of the parameters are Bjj = 62.3, 
C = 82.66 and A 3 = 2.3917. The latter corresponds to 
a value of 70 — 18.86°. The energy difference between 
the minimum of the harmonic oscillator and the values 
at one of the borders of the A7 = 60° wedge is roughly 

_g_ C (7-7b) 2 „ 0T4 

which, for a reasonable value of /?o (e.g. values smaller 
than 0.3), is higher than the excitation energy of the 
higher lying states shown in Fig. 4: E s + ~ l.bMeV. 
In order to calculate the value in Eq. (44) we used the 
Bohr-Mottelson prescription [40] for the mass coefficient, 
B m = j-^AM n Rq , where A = 2, M N is the nucleon 
mass and R = tqA 1 / 3 is the nuclear radius. Therefore 
we obtain 

h 2 4ttH 2 c 2 120.72MeV 

2BZ ~ 3M N c 2 r 2 A 5 / 3 A^~ 3 ' ( ' 

This expression, combined with the proper mass number, 
has been used in Eq. (44) . The evaluation of the depth of 
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FIG. 7: Values of 70 from irrotational fit to various osmium 
isotopes (line with circles, lower panel) and corresponding val- 
ues of the parameter C (middle panel) and Bo (upper panel) . 
Predictions obtained in [41] with a different method are also 
shown (dashed line, lower panel). 



the potential well gives an indication that the wavefunc- 
tions of the states that we have fitted are all well confined 
inside the well. In fact the true shape of the potential 
is thought to depend on 7 periodically (terms like cos 37 
are the most important), therefore the harmonic oscilla- 
tor shape can only represent a convenient approximation 
which is valid around the minimum. 

We show in Fig. 7 the values of 70 (black dots) ob- 
tained from the irrotational fit for the whole chain of os- 
mium isotopes that we have analyzed together with the 
corresponding values of the parameter C (red squares, 
middle part) and Bp (diamonds, upper part). We notice 
that, while the value of 70 decreases from the left to the 
middle and then increases, the value of the parameter C 
does not show any regular behaviour with the mass num- 
ber. When one considers the extremes of the mass chain, 
the rotor becomes softer. By applying relation (44) to 
the softest case ( 172 Os) one can check that the harmonic 
approximation is still valid. 

In Fig. 7, we also display (dashed line with blue trian- 
gles) predicted effective values of 7 using K shape invari- 
ants [41]. Despite the different technique used to obtain 
those values (which are extracted using quadrupole ma- 
trix elements and B(E2) values) they agree fairly well 
with our results. For these values, the error bars result 
from errors on the measurement of the transition rates. 
In our case an estimate of the theoretical error is quite 
difficult: in principle we can repeat the fitting procedure 
many times and we can extract a distribution of values 
of parameters. To understand if our procedure is reliable 
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FIG. 8: (Color online). Values of 70 (black dots) and of the 
parameter C (red squares) extracted from 20 independent fits 
to 178 Os. From these data one obtains an estimate of the the- 
oretical error associated with the fitting procedure described 
in the text. Mean values, (70), and mean values ± the cor- 
rected standard deviation, ((70) ± cr^ 70 )), are indicated with 
dashed and dotted lines for both quantities, respectively. 



we have repeated the fit on 178 0s 20 times. We have col- 
lected the various sets of parameters and we present in 
Fig. 8 the extracted values of 70 . It is therefore possible 
to obtain an estimate of the distribution of 70 values and 
to give mean values and corrected standard deviations as 
a simple estimate for the theoretical error. 

Hartree-Fock-Bogoliubov calculations have been car- 
ried out (Ref. [29]) in the Os-Pt region with results 
pointing toward a prolate structure with 70 very close 
to 0° and with a flat potential in 7 (for the Os nuclei 
with A = 186— 192). In such self-consistent calculations 
the potential energy surface is calculated variationally, 
starting from a given effective interaction. In most cases 
no dynamical calculations are done on top of that and 
energies are often associated with the exact value of the 
minimum of this surface. Likewise, a few excited states 
can also be predicted. Our approach starts from a given 
mathematical expression for the potential V{(3, 7), which 
is then used to solve the collective model and can pre- 
dict full energy spectra. However, the particular choice 
of V((3, 7) made at present may still deviate from actual, 
more realistic, potential energy surfaces. Besides, indi- 
cations exist for ground state hexadecapole deformation 
in the Pt-Os region, as suggested by strong E4 coupling 
arising in (a, a') and (p,p') reactions. 

Very recently a description of Os isotopes (among oth- 
ers) has been proposed within the IBM-1 [42]. The au- 
thors fit the low-energy positive-parity spectra obtaining 
a good overall description. Both their work and ours, 
which rely on models based on different ingredients, de- 
scribe equally well a good fraction of the spectral prop- 
erties of Os isotopes, although a number of differences 
clearly leave room for more detailed studies. Another 
analysis of quadrupolc moments, transition and trans- 
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FIG. 9: (Color online). Electric quadrupole moments for os- 
mium isotopes in eb. Data (red squares with error bars) are 
taken from Ref. [44] . The calculated values (black dots) have 
been obtained by adjusting an additional parameter to repro- 
duce the excitation energy of the first 2 + state. 



fer rates of (platinum and) osmium isotopes in terms of 
IBM-2 [43] shows that a reasonably good agreement with 
experimental data can be obtained with a smoothly vary- 
ing set of parameters. 

In our present study to construct solutions to the Bohr 
hamiltonian for soft triaxial nuclei, the model parameters 
used determine (i) the moments of inertia (through the 
parameter A3, or Ai for the more general fit), (ii) the stiff- 
ness parameter C for the 7— oscillatory motion and (iii) 
Bd and Ad determining the Davidson potential. There 
appear rather large variations in particular in Bp and C; 
however, the deduced value of the quantity 70 (see fig. 7 
- lower part) still exhibits a very smooth variation and 
agrees very well with results using a totally different ap- 
proach [41]. We stress that the above variations do not 
at all contradict the fact that smooth variations in the 
parameters for the IBM-2 hamiltonian result in describ- 
ing the observed data in the Os nuclei, because totally 
different models are used. 

In order to test even further our model, we present 
in fig. 9 a comparison of calculated and experimental 
static quadrupole moments of the first excited 2 + state. 
The definition of the quadrupole moment adapted to the 
present case is: 



16tt 
~5~ 



22 2 
2 0-2 



3ZRq 

47T 



(2+|K*||2+; 



(46) 



where the reduced matrix element, once we exploit the 
trigonometrical simplification for the 7 part and we in- 
sert the values of the Clebsch-Gordan coefficients, is cal- 
culated to be 
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[cos7o(a 2 — a ) + 2a a sin 70 



(47) 
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The last expression still contains the matrix element of 
[3 that depends on the Davidson ground state wavefunc- 
tion, which in turn contains the parameter Ajj . This pa- 
rameter can be fixed by requiring that eu(2 + ) — en(0 + ) 
matches with the absolute energy difference in MeV. The 
results collected in fig. 9, are correct both in sign and 
trend, but they underestimate a bit the measured values 
[44]. A more complete study, including B(E2) values, is 
in progress. 

IX. CONCLUSIONS 

In the present paper, we have presented a method in 
order to solve the Bohr-Mottelson collective model for a 
general soft triaxial nucleus. By choosing a potential of 
the form V{P,i) = Vl(J3) + V 2 (-f)/(3 2 , the Bohr Hamil- 
tonian separates exactly in the (3 and 7 variables. Using 
a displaced harmonic oscillator potential in the 7 direc- 
tion, representing the softness in the 7 variable, allows us 
to solve approximately the 7-angular equation. In solv- 
ing that part of the problem, we discuss two possibilities 
to further separate the 7 variable from the projections 
of the angular momentum. One approximation neglects 
fluctuations of the moments of inertia in the 7 direction, 
but keeps the softness. Under this assumptions, the three 
components of the moment of inertia are related and re- 
sult in one parameter. The other approximation allows 
a fitting of the three components of the moment of in- 
ertia independently. This is in line with the observation 
that the experimental data do not agree so well with the 
assumption of irrotational flow. We then present an al- 
gebraic method in order to find approximate solutions 
to the set of coupled equations that result for the mo- 
tion in the 7 variable and also discuss at some length the 
state labeling problem in order to characterize the various 
collective bands that result. We subsequently study the 
equation for the (3 degree of freedom, which can be solved 
exactly for a number of interesting potentials. We thus 
consider the Coulomb/Kratzer potential, the Davidson 
potential (which is an extension of the harmonic oscil- 
lator potential in (3) as well as the infinite square-well 
potential in our calculations. In each case, we are able to 
derive analytic expressions, which are then used to deter- 
mine the full energy spectra. These energy spectra, that 
describe oscillatory behavior in both the [3 and 7 vari- 
ables, as well as the rotational motion, in general depends 
on 6 parameters: three that characterize the V((3,j) po- 
tential and three that determine the moment of inertia. 
Using appropriate scaling in the energy spectrum only 5 
remain to be determined. Going back to irrotational mo- 
tion, the three components of the moment of inertia are 
related and this reduces the full set to just 3 parameters. 
Due to the number of parameters and the very involved 
energy eigenvalue expression, we have used a numerical 



method based on a random walk procedure in order to 
obtain the optimal fitting parameter set. We have used 
this method to study the Osmium isotopes in the interval 
172 < A < 192. One of the results is the fact that the 
equilibrium 70 value shows a rather smooth dependence 
on A, but the stiffness of the potential in the 7 variable 
indicates rigid cases for intermediate masses and softer 
cases for isotopes sitting at the extremes of the considered 
mass chain. To elucidate this point it would be very help- 
ful to search for solutions of Eq. (5) with more involved 
periodic potentials [45] . Still, the calculated values of 70 
are in rather good agreement with an independent ap- 
proach to extract an effective 7 value by Werner et al. 
We have in mind to study the rare-earth region in a more 
systematic way using the methods presented here. 

Recently other works have appeared that treat simi- 
lar problems, yet with different spirit and aims. Caprio 
[16] analyzes the numerical diagonalization of a (3— soft, 
7— stabilized problem, evidencing how the approximate 
separation of variables of the X(5) model may be ques- 
tioned on the basis of a strong (3 — 7 coupling. Despite 
the fact that our approach is complementary (rather than 
solving exactly the numerical problem, we postulate an- 
alytically solvable potentials and we try to see if they 
can be profitably used), we reach similar conclusions. 
For example we also find that a considerable degree of 
dynamical 7— softness is needed (for realistic cases of 7 
stiffness) to account for the energies of the 7— bands' lev- 
els. Furthermore our method allows to look for the op- 
timal value of the actual moments of inertia. A point of 
difference is, instead, the pattern for the staggering: in 
the examples below we rather find a scheme of the type 
(2 + , 3 + ), (4 + , 5 + ), ... similar to that of the rigid triaxial 
rotor. Another very important series of paper by Rowe 
and collaborators [13, 17, 18] furnishes a rapidly converg- 
ing method for exact numerical treatment of the problem 
by using a basis defined in terms of "deformed" wavefunc- 
tions in [3 and five-dimensional spherical harmonics. 

The authors are most grateful to J.Wood for extensive 
discussions and L.F. wishes to thank A.Vitturi for pro- 
found comments. Financial support from the "FWO 
Vlaanderen" (L.F. and K.H.) and the University of Ghent 
(S.D.B. and K.H.), that made this research possible, is 
acknowledged. 



APPENDIX: Matrix elements 

List of non-null matrix elements of operators as dis- 
cussed in Eq. (11) of Sect. II. 

{LMK I L% I LMK') = K 2 if K = K' (48) 
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f L(L + 1) - K 2 



(LMK | L\ | LMK') = 



(LMK | L\ | LMK') 



(LMK | L 2 . | LMK') 



if A ' = K 
if A"' = K - 2 
if if' = A + 2 



(49) 



where 



and 



(LMK | L\ | LMK') = 



{ L{L+l) - K 2 

2 ; 

(LMK | L\ | IMif') 
4 

(LMK | j> | LMK 1 ) 



if K' = K 
if A' = A' — 2 
if A"' = K + 2 



(LMK | | LM(K - 2)) = ^/(L + A)(L + A - 1)(L - A" + 2)(L - A' + 1) 



(LMK | l£ | LM(K + 2)) = ^/ (L + A + 2)(L + A + l)(L - K - l)(L - A) 

I 



(50) 
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